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Abstract 

We study phases and fractal structures of three-dimensional simplicial quantum gravity 
by the Monte-Carlo method. After measuring the surface area distribution (SAD) which is the 
three-dimensional analog of the loop length distribution (LLD) in two-dimensional quantum 
gravity, we classify the fractal structures into three types: (i) in the hot (strong coupling) phase, 
strong gravity makes the space-time one crumpled mother universe with small fluctuating 
branches around it. This is a crumpled phase with a large Hausdorff dimension dn — 5. The 
topologies of cross-sections are extremely complicated, (ii) at the critical point, we observe 
that the space-time is a fractal-like manifold which has one mother universe with small and 
middle size branches around it. The Hausdorff dimension is du ~ 4. We observe some scaling 
behaviors for the cross-sections of the manifold. This manifold resembles the fractal surface 
observed in two-dimensional quantum gravity, (iii) in the cold (weak coupling) phase, the 
mother universe disappears completely and the space-time seems to be the branched-polymer 
with a small Hausdorff dimension dn — 2. Almost all of the cross-sections have the spherical 
topology S*^ in the cold phase. 
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1 Introduction 



In recent years remarkable progress has been made to quantize the field theories with 
general covariance . Especially, in two-dimensional quantum gravity it is very encouraging 
that some analytical formulations (Liouville field theory [Q, matrix models^, topological 
field theories0 [||) and numerical calculations by the dynamical triangulations (DT) [0 
have shown the good agreement. It is found that the universal structure which leads to the 
continuum limit seems to be its fractal structure]^. It seems very important to investigate 
the fractal behavior of the space-time manifold in any dimension. Below, we will briefiy 
review this fractal structure in two dimension as the starting point. 

On the other hand, in the three-dimensional quantum gravity which is our primary 
concern of this paper, as well as in the four-dimensional case, we have no analytical mod- 
els which predict the scaling or fractal properties near the critical point. By numerical 
simulations, a few groups have studied the three-dimensional Euclidean quantum grav- 
ity based on dynamical triangulations and reported the existence of the first-order phase 
transition [0 |jl4|. Numerical calculations show that the entropy is exponentially bounded 



and the thermodynamic limit exists in three-dimensional simplicial gravity [T^[|l6]. The 



starting point of this model is based on the ansatz that the partition function describing 
the fiuctuations of a continuum geometry can be approximated by a weighted sum over all 
simplicial manifolds (triangulations) T which is composed of equilateral tetrahedra: 

T 

In this work the topology of the simplicial lattice is restricted to the three-sphere S^. The 
action is taken to be of the form 

So{T) = -KoNo + KsN^, (2) 

where iVj is the number of z-simplices in a simplicial lattice T, and kq, are dimensionless 
gravitational and cosmological constants on the lattice. This action (0) corresponds to a 
discretized version of the continuum action of three-dimensional Einstein gravity: 

^eh[^7] = / d^x^R + \j d'x^ , (3) 

where G is the Newton constant and A is the cosmological constant. 

To control the volume fiuctuations we add an additional term to the action of the form: 

5S = ^{Ns-Vf , (4) 

where V is the target size of the simplicial lattices and 7 is an appropriate parameter. 



In ref.[|T6| it is argued that lattices with A^3 ~ are distributed according to the correct 
Boltzmann weight up to correction terms of order 0{^^). 

In summary, our starting point is the following partition function: 

Z = Ee-^("^ (5) 

T 

where the total action S{T) on the three-dimensional simplicial lattice is of the form 

S{T) = -KoNo + K3Ns + ^{Ns-Vf . (6) 
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Now, we would like to explain the method to detect the fractal structure, taking two- 
dimensional quantum gravity as an example in some detail. Let us consider a triangulated 
surface E which is composed of equilateral triangles. One of the points to understand the 
continuum limit of the DT in two dimension is its fractal structureQ, characterized by the 
loop length distribution (LLD) defined as follows: 

1. Let us consider a 0^-graph IJ dual to the triangulated surface E. 

2. Pick up a single point P in (a triangle in S). 

3. Find all the points which has a geodesic distance (= minimum number of steps of two 
points in the dual graph) D from P. 

4. The set of links (also links in E) connecting two points of distance D — 1 and D 
consists of some closed loops in S. 

5. We define the loop length distribution (LLD) function by counting the number of these 
closed loops which have a geodesic distance D and a loop length L, and write g{L, D). 

As shown in ref.|p], the function g{L,D) is universal and explicitly calculated as 

g{L, D)xD^ = -^{x-'/' + ^x-^/^ + ^x^/^)e-^ (7) 
7\/'K / 3 

with X = L/D^ as a scaling parameter. Numerically, this fractal behavior was established 
by N.Tsuda and T.Yukawa p||Tl. 

In the next section, we will explain our data and the interpretations of the measurements 
will be described, mainly on the three-dimensional analog of the fractal structure described 
above. 



2 Numerical simulations 

To discuss the phase structure of simplicial gravity, the integrated scalar curvature per 
volume has conventionally been used. The lattice analog is: 

/ d^x^/gR (r-n(l)) Nn 



d-'x^ 



where c = 2tt/ arccos(l/3) = 5.104 ... is the number of tetrahedra around a link which gives 
the zero scalar curvature for fiat space. Here n{l) denotes the number of tetrahedra sharing 
the link /, called the coordination number. A conventional indicator of phase transitions is 
the specific heat, which in this case corresponds to the fluctuation in the number of vertices 
No: 

XN, = l^m) - {Nof) , (9) 

where averages are taken over ensembles of configurations with a fixed volume N^ = V. In 
fig.0, we show the expectation value < N0/N3 > and xnq as a function of the gravitational 
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coupling constant kq. We cannot observe the hysteresis curve, in contrast to the calculation 
in ref. [T^. We observe a phase transition near = 4.0, but with this data above, we cannot 



determine the order of the transition unless performing the finite-size scaling analysis. 
Here we should make a remark about the thermalization in our study which will explain 



the difference from the paper |]T3|. For the error estimate we used the jackknife method [pO[]. 
We examined the magnitude of error as a function of the bin size and adopted the bin size 
where the jackknife error had leveled off. This bin size corresponds to about 500 sweeps 
where one sweep is defined as the total number of accepted moves divided by the number 
of tetrahedra N^. Comparing to the typical interval of about 100 sweeps in two dimension, 
the relaxation time of the Markov chain in three dimension is much longer. We need to 
perform computer runs with the long interval larger than 500 sweeps to accumulate the 
statistically uncorrelated data. We have taken the interval 600 sweeps to calculate our 
data. The reason why the hysteresis has disappeared in fig.|I| is that we have thermalized 
our data enough to avoid the autocorrelation. In ref.[0, 10^ sweeps were taken to get 
each point in the hysteresis, but we have taken 600 x 120 = 7.2 x 10^ sweeps to get each 
point in the fig.|l|. No hysteresis is often observed in a finite system, though there exists 
a first-order transition in the system [p!7|. We cannot determine the order of the transition 



without performing the finite-size scaling analysis flS 



In general, if a statistical system has a second-order phase transition, some scaling 
invariant properties, or, fractal structures are expected to appear at the critical point. 
In the case of two-dimensional gravity, the LLD characterizes the fractal property of the 
surface [H as described in the last section. By the analogy to two-dimensional case, we 
introduce the surface area distribution (SAD) in three-dimensional case: 

1. Let us consider a connected 0^-graph T dual to a three-dimensional simplicial lattice 
T with a spherical topology 5''^. 

2. Pick up a single point P in T (a tetrahedron in T). 

3. Find all the points which has a geodesic distance (= minimum number of steps of two 
points in the dual graph) D from P. 

4. The set of links (triangles in T) connecting two points of distance D — 1 and D consist 
of some closed surfaces in T. 

5. We define the distribution function of these closed surfaces which have a distance D 
and area A as surface area distribution (SAD), and denote by f{A,D). 

Theoretically, there exist no analytical models which predict the scaling behavior of this 
SAD function f{A,D). We have measured the SAD functions in the case of the lattice 
size V = 10000. At each distance D, we compute the number of branches (baby universes) 
iVb(-D), the volume within the distance V{D) , and moreover the large area (= the cross- 
section of the mother universe) y4iarge(-D). In general, we expect these observables to scale 
as, 

V{D) ~ D'^^, (10) 
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If the manifold (space-time) has a mother body (huge universe) which dominates the fractal 
structure, we expect (in = c?m + 1- 



In fig.0, we show the data in the strong coupling limit kq = 0. There exists a huge mother 
universe. The behavior of N\^{D) in fig.H(a) is similar to that of the crumpled surface in 
two-dimensional quantum gravity From the linear region D ~ 8 in fig.|^(b), we can 
read V{D) ~ D^, or the Hausdorff dimension du ~ 5. In fig.0(c), we show the SAD with 
X = A/D^ as a scaling variable, because this variable gives the most proper fitting for the 
distribution of the large area (= cross-section of the mother universe). Fig.0(d) also shows 
that the large area should scale as yliarge(-D) ~ D^, consistent with the naive expectation 
(in = c^M + 1- In fig-@(c), we can observe no contributions from the universes of middle size. 
These facts give us the picture that the space-time manifold has one extremely crumpled 
huge mother universe which is not expanding widely and many small baby universes of 
low density are fluctuating around it. The SAD function for the mother universe which 
dominates the fractal structure of the manifold seems to indicate a scaling behavior with a 
scaling variable x = A/D^. 



Next, the data near the critical point kq = 4.0 are shown in fig.|^. The behavior of 
A^b(-D) in fig.^(a) is very similar to that of the fractal surface in two-dimensional quantum 
gravity||ll|. This is one of the reasons why we call the manifold at the critical point fractal- 
like. We can read dh — 3 from fig.^(a). From the linear region in fig.|^(b), we can read 
the Hausdorff dimension c/h — 4. Our data show the existence of the mother universe, 
which extends more widely than that of the hot phase. Fig.|^(d) shows that the large area 
(= cross-section of the mother universe) scales as Ai^rge ~ D^, giving the same Hausdorff 
dimension (in — 4. So we choose x = A/D^ as the scaling variable, but the small and 
middle areas (= rest parts after subtracting the contribution of the mother) do not scale 
as well by this variable in fig.0(c). Then, we divide the SAD function into two parts, the 
contribution of the large area /large (^; D) and rest of areas (hereafter we abbreviate small 
and middle areas to small areas) /smaii(^, D) : 

f{A, D) = /iarge(A, D) + /,^all(A, D) . (11) 

The divided data in fig.|^(b) tells us that x' = A/D'^ is a good scaling variable for the 
dominant part of /small (^, D). 

Let us discuss the fractal(-like) structure at the critical point in some detail. In fig.H(a), 
we can assume the large area distribution function as, 

/iarge(^, D) = const. X x'^e""" . (12) 

Here a and b are some constants. The fittings in the fig.^(a) show that we can safely assume 
a > —1, so the integration of /iarge(^5 D) given in eq.(|12|) is convergent. In other words, 
the following quantities are finite: 

(A")iarge = lim / rfA A"/large(^, D) {u > 0) . (13) 

Here e is the cutoff for the area on the simphcial lattice. It is easy to show that 

(AO)iargc oc r{a+l)xD\ 

(Al)large OC F {a + 2) X . (14) 
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From eg . (p^ , the large area part (mother universe) is independent of the cutoff e, so this 
part seems to be the universal quantity. 

On the other hand, from fig.^(b) the small area distribution function is parametrized 

as, 



fsmAA D) = const. X x'-"'e-'''^' . (15) 

Here a' and h' are positive constants. From the gradient of the dominant line in the fig.^(b), 
we can read the value of the exponent a' = 2/5. The good fitting function (after rescaling 
b'x' — s> x') is the following form: 

U^n{A, D) = const, x x'-^/^^-.' _ ^^g) 
Then, we can calculate the quantity: 

/oo 
dA/small(^, D) 

OC const. X D^e"^/^ + const, x Dt'^'"^ + const, x , (17) 

/oo 
dAAU,,,{A, D) 

oc const. X D^e"^/^ + const, x D"^ . (18) 
Moreover, the average area of the small part behaves as 

small ~ t . (19) 

This relation (|1^) is quite similar to the two-dimensional fractal surface where the small 
loop has a constant length with the size of the cutoff length]^. The small baby universe is 
not a universal quantity but also a local quantum fluctuation. If we identify 

(A°)small = iVb(D) ~ Z}'^\ (20) 

then we can decide (it = 3 by comparing eq. (p!7|) and eq.(^). This value of = 3 accounts 
very well for the behavior of Ni,{D) in fig.^(a). Moreover, eq.(|TB|) and eq.(|Il|) give the same 
Hausdorff dimension = 4, accounting for the behavior of V{D) in fig.^(b). 

Here we mention the observation by Agishtein and Migdal]^. It was found in ref.0 
that the Hausdorff dimension grows logarithmically 

c^H ~ ci + C2 log D (ci, C2 are constants) . (21) 

They argued that the logarithmic growth of Hausdorff dimension might be a universal 
law of Nature. But, in our interpretation, this behavior (^) expresses the non-universal 
(cutoff-dependent) property which originates in the small part (see eq. (P^)). 

Now we can explain the notion of the fractal-like manifold. At the critical point, the 
mother, middle and small universes contribute to the SAD function as well as to the LLD 
function of the fractal surface in two-dimensional quantum gravity [§]. But the distinct 
point is the observation that we have two kinds of scaling variables, x = A/D^ for the 
mother and x' = A/D"^ for the baby universes in three-dimensional quantum gravity. This 
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is why we call such a manifold not fractal but fractal-like. It is a remarkable property of 
three-dimensional fractal-like manifold that both the large and small areas give the same 
Hausdorff dimension c/h — 4, though the large area gives cin = 3 but the dominant small 
areas give du = 4 in two-dimensional fractal manifold]^. 

Finally, we show the data in the weak coupling phase (cold phase) = 5.0 in fig.^ 
In this phase, the space-time manifold consists of widely expanding branches and has no 
mother universe. Fig.|^(a) shows N];,{D) ~ D. The behavior of N)^{D) is quite different from 



the other cases, but similar to that of the branched-polymer in two-dimensional gravity] 11]. 
From fig.^(b), we can easily read that V{D) ~ D^, or c/h — 2. Because the mother body 
does not exist in this phase, we should choose a scaling variable x = A/D for the dominant 
small fluctuations in fig.|](c). 



3 Summary and discussion 

There exist two phases, hot and cold phase, in three-dimensional simplicial quantum 
gravity, as previously known |T2|[p^. They argued the first-order transition between these 



two phases because of the existence of a large hysteresis |]T3|. In our study we observed a 
rapid transition near kq = 4.0 but no hysteresis, so we cannot determine the order of the 
transition without the finite-size scaling analysis. 

By measuring the SAD functions using the Monte-Carlo technique, we could characterize 
the phases of three-dimensional simplicial gravity into three types: 

1. The hot phase (strong coupling phase). 

This phase is a typical crumpled phase. Strong gravity makes the space-time a 
crumpled manifold with small fluctuations around it. The fact that the area of the 
cross-section of the mother universe scales as y4iarge(-D) ~ Z)^ determines the Hausdorff 
dimension to be (in — 5. The topology of the cross-sections for the mother universe 
is highly complicated, though those of baby universes are almost spherical topologies. 
Because the contribution of the baby universes to the SAD function is quite different 
from that of the mother in contrast to the LLD in two dimension, we do not know 
how one associates a fractal structure with this case. 

2. The critical point. 

At this point, the space-time is a fractal-like manifold. By the term fractal-like, we 
means that this manifold seems to have two different scaling variables, x = A/D^ for 
the mother universe, and x' = A/D"^ for the baby universes. Both the long and short 
distance regions give the same Hausdorff dimension (in — 4. The fractal behavior of 
the long distance part seems to be universal or cutoff- independent. The topologies of 
the cross-sections are less complicated than those in the hot phase. 

3. The cold phase (weak coupling phase). 

In this phase, the mother universe, which exists in the hot phase and at the critical 
point, disappears completely. Weak gravity cannot make any large body in the space- 
time. The fractal structure is determined by the dominant small fluctuations with 
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a scaling parameter x = A/D. Almost all of the cross-sections have the spherical 
topology S"^. The space-time is a branched-polymer with a small Hausdorff dimension 
dn ~ 2. We could not find any universal (cutoff-independent) property in this phase. 

We should discuss the properties at the critical point in more detail. Because three- 
dimensional Einstein-Hilbert action (|^) includes a dimensionfuU coupling constant G, this 
theory is perturbatively unrenormalizable as well as four-dimensional Einstein gravity. In 
other words, we cannot expect the scale invariance over all order of scales in such a theory. 
But our data shows that two different types of scaling regions seem to exist near the 
critical point, short distance region (small baby universe) and long distance region (mother 
universe). Short distance region has its own scaling with a scaling variable x' = A/D^. 
But this region seems to be cutoff-dependent, so we cannot expect the universal structure 
in this region. On the other hand, long distance region also has its own scaling with a 
scaling variable x = A/D^. As the scaling of this region seems to be cutoff-independent, 
this fractal-like structure may be a universal one which leads to the continuum limit as 
well as the fractal structure does in two-dimensional lattice gravity. This point is very 
important. To confirm our conjecture, we added the irrelevant discretized i?^-term to the 
Einstein-Hilbert action (^ on the simplicial lattice. Though the SAD function of the small 
part changed its form, that of the mother part did not change its form when the coupling 
to i?^-term is small |jl9|. But these facts do not establish the universality of the fractal- like 



SAD function, so we need further studies based on the method in constructive field theory 
and interesting results are expected. 
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Figure 1: Cooling and heating processes with lattice size V = 10000. (a)Scalar curvature per 
volume .vs. kq for each process. No hysteresis is observed. (b)Node susceptibility .vs. kq 
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Figure 2: The data in the hot phase (strong couphng hmit kq = 0). (a)Number of branches at 
each distance D. (b)Volume within the geodesic distance D. (c)The SAD is shown with a scahng 
parameter x = A/D'^. (d) The expectation value of the large area. 
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Figure 3: The data at the critical point {hq = 4.0). (a)Number of branches at each distance 
D. (b)Volume within the geodesic distance D. (c)The SAD is shown with the scahng parameter 
X = AjD^ . (d) The expectation value of the large area. 
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Figure 4: The divided SAD function at the critical point (kq = 4.0). (a)The SAD for the large 
part measured and some fittings. (b)The SAD for the small part measured and some fittings. 
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Figure 5: The data in the cold phase (kq = 5.0). (a)Number of branches at each distance D. 
(b)Volume within the geodesic distance D. (c)The SAD is shown with the scahng parameter 
X = A/D. 
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